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Abstract 

We prove the existence of quadrature formulas exact for integrating high degree polynomials 
with respect to Jacobi weights based on scattered data on the unit interval. We also obtain a 
characterization of local Besov spaces using the coefficients of a tight frame expansion. 



1 Introduction 

It is well known that a major drawback of polynomial approximation is that polynomials cannot be 
localized; a polynomial of degree n is completely determined by its values at n + 1 points on an interval, 
howsoever small. Another example is the following. If P* is a best polynomial approximation of degree 
at most n to the function f(x) = \x\ on [—1, 1], then there are at least n + 2 points yj on [—1, 1], where 
\f(Vj) ~ Pn(Vj)\ — cn _1 , where c is a positive constant independent of n and the points yj. Moreover, as 
n — > oo, these points become dense on [—1, 1]. Thus, even though the function is piecewise analytic, the 
only singularity of the function, namely, x = 0, affects the quality of best approximation on the entire 
interval. 

One of the interesting problems in the theory of polynomial approximation is thus to construct local- 
ized polynomial approximations. There are two flavors to this problem. In the local information problem, 
one has information about the target function only on a part of the interval [—1,1], and wants to con- 
struct the approximation also on a neighborhood of this part, keeping the growth of the approximation 
in control on the rest of the interval. In the local adaptivity problem, the data is available on the whole 
interval; indeed, it might even be in the form of the globally defined Fourier coefficients f(p; k) (defined in 
(|2.4p ) of the target function / for a suitable measure /z. Nevertheless, we wish to find a single polynomial 
that adapts its behavior on different parts of the interval according to the smoothness of the target func- 
tion on these parts. This is a different problem from piecewise polynomial (or spline) approximation and 
adaptive approximation. Unlike in spline approximation, we want to obtain a single polynomial, which 
yields the benefits of the superior degree of approximation provided by polynomials on intervals where 
the function is smooth or even analytic. Unlike in adaptive approximation, we do not wish first to find 
the location of the singularities of the target function; the determination of their locations should be a 
byproduct of successive polynomial approximations. The construction of these polynomials may depend 
upon finitely many of the Fourier coefficients, or on the values of / at certain points on the interval. In 
most modern applications, one does not have a control on where to choose these points. Such problems 
are called scattered data problems. Thus, one of the important problems in this theory is to develop 
quadrature formulas similar to the Gauss-Jacobi quadarture formula, except that one should not require 
that the values of / be at the zeros of some orthogonal polynomial, but may be at any set of points 
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instead. In return, we cannot expect the formula based on m points to be exact for polynomials of degree 
2m — 1 but only for polynomials of degree proportional to m. 

To summarize, we wish to construct a wavelet-like expansion X^fe a j,fe(/)V ; j,fe with the follow- 
ing properties: (1) For k = 0, l,---, V^fc is a polynomial of degree at most 2 k , independent of /. 
(2) Each coefficient a J; fc should be a finite linear combination of the Fourier coefficients of the target 
function /, or the values of / at C(2 fc ) arbitratily chosen points on the interval. (3) The difference 

Hk=o 52j a jAf)^jA x ) ~ f( x ) should be 0(E^ c2 n f00 (f)) uniformly on [-1, 1] (sec (J2ZTJ for definition, 
fi being a suitable measure), and yet, near every point x € [—1,1], it should be small according to the 
behavior of / near x. (4) It should be possible to characterize the smoothness of / near different points 
by the absolute values |»j>(/)|. (5) The following frame inequalities should be satisfied: 

EEl fl ^(/)l 2 HI/IU;2- (1.1) 

fe 3 

The system {ipj^k} will be called a frame. In the case when equality holds in (|1.1[) rather than ~, the 
system {V'j.fe} will be called a tight frame. 

This paper is a continuation of our paper [13], where we constructed certain localized kernels, and 
characterized local Besov spaces on the unit interval in terms of polynomial frames based on Jacobi poly- 
nomials. The aim of this paper is to fill in certain gaps left over in [13] . which have come to our attention 
after the publication of that paper. In Section^ we introduce the necessary notations. In Section [3] we 
develop the localized tight frames in the context of general mass distributions on [—1,1]. Essential ingre- 
dients in this theory are localized polynomial kernels and quadrature formulas. In Section 3] we develop 
quadrature formulas exact for integrating high degree polynomials with respect to the Jacobi weight, 
based on scattered data on the interval. In Section [5] we prove localization estimates on certain kernels 
based on Jacobi polynomials that are more elegant than those presented in [13j . Certain computational 
issues are addressed in Section [5] 



2 Notations and background 

For x > 0, let II^ denote the class of all (algebraic) polynomials of degree at most x. Denoting by [x\ the 
integer part of x, this is the same as the class II y . Extending the notation in this way allows us to use 
a less cumbersome notation, as well as relieves us of stating certain conditions on the degrees in almost 
every statement. If [i is a finite positive or signed Borel measure on [—1,1], the total variation measure 
of /i will be denoted by A point x E [—1, 1] is called a point of increase of \i if \fi\(I) > for every 
open interval I containing x. The set of all points of increase of /i is called the support of fi, denoted by 
supp (fx). It is clear that supp (/i) is always a closed set. The measure /i is called a mass distribution if 
it has infinitely many points of increase and |/x|([— 1, 1]) < oo. If A C [—1, 1] is a /i-measurable subset, 
and / : A — ► C, we define 



l/IU;p,- 



A 



i/p 

\f(t)\ p d\n\(t)} , ifl<p<oo, 



\fi\ - ess sup teA \f(t)\, ifp = oo. 



The class L p (fi,A) consists of all /^-measurable functions / : A — ► C for which H/H^^a < oo, with 
the convention that two functions are considered equal if they are equal |/^|-almost everywhere. The 
space C(A) denotes the class of all uniformly continuous, bounded functions on A. The symbol X p (fi, A) 
denotes L p (fi,A) if 1 < p < oo, and C(A) if p = oo. For / S L p (fi, A) and x > 0, we define the degree of 
approximation of / from by 

E mA (f) := m^Jf - P\\ KP , A . (2.1) 

In the sequel, if A = [—1,1], we will often omit its mention from the notations, when not required by 
considerations of clarity. Similarly, we will often omit the mention of the measure fj, if it is the Lebesgue 
measure. 
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It is convenient to define the Besov spaces using the degrees of approximation. Characterizations in 
terms of the Ditzian-Totik moduli of smoothness are given in [5] for some measures p. Let < p < oo, 
7 > 0, and a = {a„}^L be a sequence of real numbers. We define 

oo 

{J22 nJP K\ p } 1/p , if0<p<oo, 

n=0 (2-2) 
sup2™ 7 |a„|, if p = oo. 

n>0 

The space of sequences a for which ||a|| Pi7 < oo will be denoted by b pn . For 1 < p < oo, < p < oo, 
7 > 0, the Besov space B KP ,p, 7 := B^-p^ consists of functions / G X p (p) for which the sequence 
{En;2™, P } S b p , 7 . We define 

ll/lb, !P , P „ = ll/t;P+ll{^:2",p}llp,7- (2-3) 

Let xo G [—1, 1]. The local Besov space B ll - PiPn (xo) consists of all / S X p (p) with the property that for 
every C°° function <f> supported on an interval containing xq, f4> G B^ pp/1 . Let p be a finite positive 
measure on [—1,1]. If p is a mass distribution, one can use the Gram-Schmidt orthogonalization process 
to obtain a unique sequence of orthonormalized polynomials pk(p;x) = r yk(p)x k + ■ ■ ■ G life, Jk(p) > 0, 
k = 0, 1, • • •, such that 

l 

Pk{p;x)p {p;x)dp{x) = 0, if k^j, k,j = 0, 1,- • ■, 

•l 

and j\pl(p; x)dp(x) = 1, k = 0, 1, • • •. It is customary to define Pk(p~, x) = if k < 0. If / € ^(p), we 
may define 

f{p;k) = J J(t) Pk (p;t)dp(t), k = 0,l,---, (2.4) 

and for m = 1, 2, • • •, 

m — 1 „i 

s m (p;f,x) := f(p;k)p k (p;x) = f(t)K m (p;x,t)dp(t), (2.5) 
fe=o 

where the Christoffel-Darboux kernel K m (p;x,t) :— X^fcLo 1 Pk(p',x)pk(p',t). In particular, 
l 

K m (p;x,t) 2 dp(t) = K m (p;x,x), x G K, m=l,2, (2-6) 

l 

We find it convenient to introduce the notation 

X m {p'i x)^ 1 := K m (p; x,x), x G M, m = 1, 2, • • • . 

In the sequel, the symbols c,c\, ■ • ■ will denote generic positive constants independent of the degrees 
of polynomials involved and the target function /, but may depend upon such fixed quantities in the 
discussion as p and the smoothness parameters 7, p, etc. The symbol A\ ~ A 2 means that C\A\ < A 2 < 
c 2 Ai. 



3 Polynomial frames 

We will assume in the sequel that p is a fixed mass distribution. 

Let S > 2 be an integer, and h° : R — > [0, 00) be a compactly supported function that can be expressed 
as an S times iterated integral of a function of bounded variation on R. In view of the Poisson summation 
formula [T], 

^2h°{k/\y kx = I h°(t)exp(i\{x + 2kir)t\dt, x G [-%, vr], A > 0. 
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It is not difficult to derive from this formula by a repeated integration by parts that (cf. [15] ) 



< C y(^ o(S_1) )Amin{l,(A|a;|)- 5 }, x G [-7r,7r],A> 0, 



where V^(/i ' s_1 - ) ) is the total variation of h ^^ 1 ^ on M. The following definitions are motivated by this 
observation. 

If h : [0, oo) — > R is compactly supported, let 

<S>x(ti;h,x,t) :=22 h T ^(^^(k*). x,teR, A > 1. (3.1) 



fc=0 



If A < 1, then we set $a(m; x i t) — 0- It may seem strange to define it so. However, in most applications, 
h(t) = if t > 1. Therefore, even if <&\(fj,;h,x,t) were defined by (|3.1[) for all A > 0, we would have 
<$>\(/j,;h,x,t) = h(0)po(fi) 2 — $i(fj,;h,x ) t) for all A < 1. The statements of our results later (e.g., 
Corollary 13. 2j) will be simplified with the definition as we have made it. We note that if we choose 
h(t) = 1 for < t < 1 and h(t) = if t > 1, then for integer n > 1, the kernel ^) reduces to the 

Christoffef-Darboux kernel. In the sequel, we will assume that A > 0. 

Let S > 1 be an integer. A function h : [0, oo) — > [0, oo) will be called a multiplier mask of order S if 
each of the following conditions (|3.2[) . (13. 3p is satisfied: 

sup / |$a(a*5 h,x,t)\dfi(t) < oo, (3.2) 
a>o, xe[-i,i] •'— l 

sup |$ A (/i;/i,x,i)| < c(^)A~ s , A,(5>0. (3.3) 

If v is a (possibly signed) measure, we define the summability operators by 

o x (v,n;h,f,x) := J $ x (n;h,x,t)f(t)dv(t), A > 0, x G R. (3.4) 
If ^ = /i, then we write cr^(/z; /i, /) := ^(f, /). We note that if A > 1 then 

<r\{p,h,f,x) = I - j f(fj,;k)p k (fJ,;x), (3.5) 



fc=0 



Since h is compactly supported, c^(/i; h, /, a;) can be computed using 0(A) Fourier coefficients of /. If v 
is a discrete measure that associates the mass w z with each point z in a finite set C, then 

a\(v,pL\h,f,x) = w z ®\([i; h, x, z)f(z), 
zee 

which can be computed using the information {f(z)} z ^c- Our theory usually does not depend upon the 
choice of the set C and the specific weights w z . Therefore, we will often use the Stieltjes notation to 
denote a sum of the form J^zec w zf( z )- Thus, we define a measure v to associate the mass w z with each 
z G C, and write J fdv or a variant thereof to denote the finite sum. If 1 < p < oo, then we will say that 
v < P M if / G XP(p) implies that / G and \\f\\ v;p < c||/|| KP . 

A (possibly signed) measure v will be called a M-Z (Marcinkiewicz-Zygmund) quadrature measure 
of order n (for [i) if each of the following conditions (13. 6|) . ()3.7|) is satisfied: 

y Pd^ = J Pdfi, P G n„, (3.6) 

||P|U;l <c||P|U;i, PGn„. (3.7) 

(As usual, we tacitly think of the measure v to be a member of a sequence of measures.) 

In the rest of this section, let h : [0, oo) — > [0, oo) be a nonincreasing function, such that h(t) = 1 if 
< t < 1/2 and h(t) = if t > 1. We will assume that /i is a multiplier mask of order S 1 . The constants 
c, Ci, • • • may depend upon /i. 

The following proposition lists some immediate properties of the summability operators. 



4 



Proposition 3.1 Let m > be an integer, v be an M-Z quadrature measure of order 3m/2 — 1. Let 
l<p<oD and f G L p (jx). 

(a) For P G U m/2 , a m (v, //; h, P) = P. 

(b) We have 

\\a m (u,n;h,f)\\tt., p <c\\f\\ v . tP , m = 0,l,---. (3.8) 
Consequently, if v < v fj, and f G X p (/i), then \\a m {v, At; h, f)\\p,;p < c ll/IU;p; anc ^ 

E^, m ,p(f) < \\f ~ o- m (v, n;h, f)\\^ p < cE fl . m/2: p(f). (3.9) 

(c) // / is supported on a subinterval I of [— 1, 1], and J is an interval with I C J C [—1, 1], i/ien 

||er m (i/, A*5 /i, /)IU;oo,[-i,i]\J ^ c||/|| M; im" s ', (3.10) 

where c may depend upon I and J in addition to fi, S, and h. 

We point out an application of this proposition for the local information problem. Suppose J is a 
subinterval of [—1, 1], C is a finite subset of J, and we have to construct an approximation to / G 
based on the values {/(z)} z6 c- We may take a subinterval / of J, and a C°° function which is equal 
to 1 on / and outside J. The above proposition shows how to approximate f on I while keeping the 
growth of the approximating polynomial under control. We find an integer m such that there exists an 
M-Z quadrature measure v m of order 3m/2 on a set containing C as its subset, and v m fi. The 
existence of such measures is proved in |18l Proposition 5.1] under the assumption that fi({x}) = for 
each x G [-1,1]. Then ([33)1 shows that 

11/ - o- m (v m , fj,; h, /0)IUoo,/ < cE^ 

;m/2,oo 

while (|3.10p and (|3.8|) limit the growth of \\o- m {v m , \x; h, /</>)|| W oo,[-i,i]U- 

Proof of Proposition 13.11 If P G U m / 2 , x G M, then $ m (fi;h,x,o)P e n 3m/2 _i. Since ft(t) = 1 if 
t < 1/2, it is easy to verify using (|3.6[) that 

$ m (^; ft, a:, t)P{t)dv{t) = f $ m (n; h, x, t)P{t)d^{t) = P(x). 
-i J-i 

This proves part (a). 

In view of (|3.7[) and (|3.2p , we have 

sup / \^ m (ji',h,x,t)\d\v\(t) < c sup / \$ m (jj,;h,x,t)\d(i(t) < c. 

xE[-l,l]J-l xE[-l,l]J-l 

The estimate p.8[) is now a simple consequence of the Riesz-Thorin interpolation theorem and the fact 
that v < p (x (cf. [T51 Lemma 4.1]). 

If P G n„j/ 2 is arbitrary, then part (a) and (|3.8[) imply that 

E nm ,p{f) < \\f - <j m (v,fi;h,f)\\^ p = \\f - P -a m (v,fj,;h, f - P)\\ KP < c||/-P|| M;p . 

This proves ([3T9"]) , 

The estimate (|3.10[) is easy to deduce using (|3.3[) . □ 

Corollary 3.1 Let m > be an integer, v be an M-Z quadrature measure of order 3m— 1, and v /i. 
Then 

H-PIUp ~ \\p\U-, p , p& n m , l < P < oo. (3.ii) 

Proof. If / G L 1 ^), then J3T|) and (applied with yU in place of v) imply that 

lk!m(M; < c||o-2 m (M; K /)IU;i < c ll/IUi' 
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We note that v ~<oo [i and that a\ m {p,\ h, f) £ C([— 1, 1]). An application of (13. 8j) again with p — oo and 
/i in place of ^ shows that 

IKm(M,/)IUoo < c||(T 2 * m (^;/i,/)|| M;00 < c||/|| woo , / G 

Consequently, the Riesz-Thorin interpolation theorem yields ||02 m (/i; ^ /)IU;p — c II/IIm;p f° r every p £ 
(l,oo) and / £ L p (/i). Using this estimate with P £ H m in place of /, we obtain from part (a) of 
Proposition 13 . 1 1 that ||-P||„ :p < c||P|| M;p , 1 < p < oo. The estimates in the reverse direction follow directly 
by using (|3.8|) with P in place of / and recalling part (a) of Proposition 13.11 again. □ 

The following theorem describes the local adaptivity property of the summability operators. 

Theorem 3.1 Let 1 < p < oo, f £ X p (/j,), x a £ [-1,1], < p < oo, 7 > 0, S > max(l,7). For 
each integer n > 0, let v n be an M-Z quadrature measure of order (3/2)(2 n ) — 1, and v n -< v ji. Then 
f £ B^. iPjPi j[xo) if and only if there exists a nondegenerate interval I centered at xq such that 

\\f - <r 2 n(v n ,iJ,;h, f)\\ K pj £ b pn . (3.12) 

Proof. Let / £ B p . P:Pn (xo), and J be an interval centered at xq such that -E M; 2n )P (/</>) £ b Pi7 for every 
C°° function <fi supported on J. We take / C I\ C J to be intervals centered at xq, and a C°° function 
V> supported on J and equal to 1 on I\ . Then (|3. 10[) leads to 

H»(i/ n , W h, (1 - V)/)IU;p,/ < J, /)2-" S . 
Since ip(t) = 1 for t £ I, we conclude that 

||/ - 0-2" {Vn,W K /)IUp,J - 11^/ - CT 2" (Vn, A 4 ) ^ V'/)IU;p,I + Ik2» (^n, (1 ~ ^)/)IUp,J 

< c (/,j,/){£; ai;2 „- 1 , p (v/) + 2-" s }. 

This proves ([535]) . 

Conversely, let (j3 . 1 2[) hold, and be a C°° function supported on /. The direct theorem of approxi- 
mation theory shows that there exists R £ II 2 r> such that \\cj) — R\\oo < c(</>)2~ . Therefore, using (|3 . 8(1 
and the fact that v n < v [i, we derive that 

E^2^, P {f4>) < \\f4>~ -Ro-2"(^n,M;^)/)IU;p 

< ||(/ - m; ft 5 /))0IU:p + 11(0 - i? ) <j 2" On, m; ^ /)IUp 

< c(0){||/- ( 7 2 „K, M ;/ l ,/)|U ;p , / + 2-" s ||/t ;p }. 

In view of (|3.12[) , this implies that / € B^p pn (xo). □ 

The discrete Hardy inequality [2] Lemma 3.4, p. 27] implies that if {aj}°^ € b Pi7 then {X^fc a j}fcLo 
b Pi7 as well. Therefore, we can deduce the following corollary using (|3.9p and (|3 . 12|) . 

Corollary 3.2 With the set up as in Theorem \3.1[ we have 

00 

/ = (o> On, m; ^ /) - ct 2"-i m; ^ /)) ( 3 - 13 ) 

im'i/l convergence in the sense of X p (p). Moreover, f £ B^ p ^ pn {xo) if and only if there exists a nonde- 
generate interval I centered at xq such that 

\\azn(v n ,fj,;h,f) - a- 2 «-i(f n -i,M;fr>/)IUp,-r e b p-.f 

The frame properties of the operators / 1— * cr 2 ™ (y n , fi; h, /) — <T 2 r>-i (v n -i, Mi /) were described in [T3] , 
We will now describe a tight frame, and demonstrate its use in characterization of local Besov spaces, 
thereby achieving all the objectives listed in the introduction. 

Let git) := yj h(t) — h(2t). Then g is supported on [1/4, 1], We assume that g is also a multiplier 
mask of order S. We define 

<G«;fr,/) := o"2»(m;5:/)- 

In the following theorem, we can choose v n to be a discretely supported positive measure. In the case of 
the Jacobi weights, it may be the measure that associates with each zero £fe, c (2 n )(/■*)) k = 1, • • • , c(2") of 
Pc(2 n ) the mass A c ( 2 n) (/i; Xfc, c (2™) (m)) f° r a suitable integer constant c ( [21[ Theorem 25, p. 168]). 



G 



Theorem 3.2 Let 1 < p < oo, f £ X p (fi), x £ [-1,1], < p < oo, 7 > 0, S > max(l,7). For each 
integer n > 0, let v n be a sequence of measures with each v n satisfying 



Pdv„ 



Pen 



2 »+l_l. 



(a) We have, with convergence in the sense of X p , 



(b) In the case when p = 2. 



CO „i 

f = 12 T*([i;h,f,t)<S> 2n (fi;g,o,t)dv n (t). 



\f\\%2 = J2\\<(^h,f)\\t n . 2 . 



(3.14) 



(3.15) 



(3.16) 



n=0 



(c) Let v n be an M-Z measure of order 3(2") — 1 and v n ~< v ji. Then f £ B^p tPt y(xo) if and only if there 
exists a nondegenerate interval I centered at Xq such that {||t*(/z; h, f)\\u-p,i} G bp, 7; or equivalently, 
there exists a nondegenerate interval I centered at Xq such that {||t*(/z; h, f)\\v n ; P ,i} G b Pj7 . 

We remark that by choosing v n to be a discretely supported positive measure, the representation 
(1 3 . 1 5 [) is actually a double sum representation, similar to classical wavelet expansions. The equation 
()3.16p then appears as the analogue of the classical Parseval identity for trignometric Fourier series. 
Part (c) shows that the absolute values of the coefficients (t£(/x; h, /, i)) tg supp 0„) m the representation 
()3.15j) can be used to characterize the local Besov spaces. It should be noted, however, that the functions 
3>2™(M;g,o,i), t £ supp (f n ), are not necessarily linearly independent. 

Proof of Theorem 13.21 In view of (|3. 14|) . a straightforward calculation using the orthonormality of 
the polynomials Pk{^) shows that for i£R, 



t„ (m; h, /, t)$ 2 n (fi;g,x, t)dv n {t) 



T* n (m; h, /, i)$ 2 ™ (a*; 9, x, t)dn(t) 



(3.17) 



= 0-2* {n; h, f,x) - crZ„-i(jj,;h,f,x). 
The equation (|3.15p follows from (|3. 13|) . used with each v n = fi. Similarly, 

\K(^ h J)\\l n -2 = IK(W^i/)ll£;2> 

and (|3. 16|) follows by a straightforward calculation (cf. [TBI Theorem 3].) Next, we prove part (c), where 
stronger assumptions on v n arc made. Since g(t) = if t < 1/4, and g is a multiplier mask of order S, 
we have for any P £ iI 2 n-2, 

h, f)\U P = IK(m; K f - P)IU <4f- p\\ kp . 

Thus, ||t*(/i; h, /)|| P;P < cE^n-s p (f). Therefore, if / £ B pj - iPjPn (xo), then an argument similar to 
the proof of Theorem 13.11 using the assumption that g is a multiplier mask of order S, implies that 
{\\T*((j,;h, f)\\p, p ,i} G b Pi7 for some interval I containing x . Conversely, let {||r*(^; h, G b P:7 , 

J C / be a proper subinterval centered at Xq, and (f> be an arbitrary C°° function supported on J. 
Denoting by X the characteristic function of /, we obtain in view of (13. 2j) with g in place of h that 



T n (m; K f, 0*2" (m; 9, °, t)dn(t) 



H;p,J 



< 



T*(p; h, /, t)X(t)$2« (/*; 5, °, *)<*A*(*) 



< 



c||<(m;^/)IU;p, 



(3.18) 



If a; € J, then 



/i, /, i)3 2 «(w 5, s, t)|d M (*) < c(/, J)2-" s ||<( At ; h, f)\\,, p < c(I, J)2-« a ||/|| M;p . 



)-nS| 
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The above two estimates and (|3.17[) show that 

\\a^ n (fJ,;hJ) - (T2n-i(^,h,f)\\^j = 



i~n (m; h, /, t)$ 2 » {p; g, °, t)dfj,(t) 



< c 1 ||<( / ,;/ l ,/)|| wp , / + C (/,J)2-" s ||/|| Ai;p . 
Thus, \\a^ n (p,;h,f) - <7£ n _i(/i; h, f)\\ KP ,j G b P;7 , and Corollary [321 implies that / G B KPyf>n (x ). 



□ 



4 Quadrature formulas 

In this section, we restrict ourselves to the case of Jacobi weights, and prove the existence of quadrature 
formulas at arbitrarily chosen points on the interval, exact for integrating high degree polynomials with 
respect to the Jacobi weights. 

For a, (3 > —1, let w a ^(x) := (1 — x) a (l + x) 13 , x G (—1,1), and w a ^{x) = otherwise. We define 
dp( a, P\x) = w a ^(x)dx, and denote the corresponding degree k orthonormalized polynomial Pkip^ 01 '^) 

with positive leading coefficient by pi™ , (k = 0, 1, ■ • ■). In particular, for integers k,j — 0, 1, • • •, 

l i ^(*)p i ^(«)«, ai ,(*)d B ={ j; 2 e ZL 

In the case a = f3 = —1/2, one obtains the Chebyshev polynomials. For 9 G [0, tt], let Tk(cos9) := 
cos We have for x G [-1,1], p { ~ 1/2 '~ 1/2) (x) = [l/^)T (x), and for k = 1, 2, • • • , p { ~ 1/2 '~ 1/2) (x) = 
(^2/^)T k (x). 

Let M > 4 be an integer, = 9q < Q\ < ■ ■ ■ < 9k < 9k+i < • • • < 9m < 9m+i = tt be arbitrary points, 
and Zk — cos 9k, k = 0, • • • , M+ 1. We will tacitly assume that the set C = {zk} is one of the members of 
a nested sequence of sets of points on [—1,1], whose union in dense in [—1,1]. Thus, all constants will be 
independent of M and the points Zk- The mesh norm (also known as fill distance) of the set is defined 
by 

8c ■= max min \9 - 9 k \ = - max \9 k +i - k \. 

0e[O,7r] l<k<M 2 0<fe<M 

The separation radius of C is defined by qc := (1/2) min <fe<M |#fc+i — 6k\- The set C is called p-uniform 
if 8c < pqc, and uniform if it is p-uniform for some p. Again, it is to be understood that the value 
of p is the same for the whole implicitly understood sequence of sets C. If 8c < 7r/4, and m > 1 

is the integer part of 7r<5 r T 1 /4, then each interval of the form [ — — , - - — 1 contains at least 

4m 4m 

one 9k- We choose one 9k j from each such interval to obtain a subset C = {cos 6^} C C such that 
Qc > 7r/(4m) and 8c < 5c < n/m < 48c- Moreover, {9 kj } C [c/m, n — c/m]. Thus, in the sequel, we 
may assume that the sets C are all 4-uniform, and that there exists a constant c with the property that 
{cos -1 z : z G C} C [c/m, tt — c/m]. Clearly, for a uniform set C, qc ~ <5c ~ 

Theorem 4.1 Let a,/3> —I, p > 1. There exists a constant a > depending only on a,{3,p, with 
the following property. Let C C (—1,1) &e a finite p-uniform set, and m > 1 be an integer such that 
7r/(2m) < (5c < n/m. Then there exist positive numbers w z , z G C, with the following properties: 



(4.1) 



^w z p{z)^ [ P(t)w a ,p(t)dt, Pen am , 

zee J- 1 

w z ~\ m {^ a > 0) ;z), zeC, (4.2) 

$3«;,|P(*)| ~ f \P(t)\w a ,p(t)dt, PeIL am . (4.3) 
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The proof of this theorem follows a familiar theme, introduced in [3 [14]. We first recall an abstract 
quadrature formula in the setting of general finite dimensional spaces (cf. [HJ Theorem 3.2.1]). 

Proposition 4.1 Let (X, \\ ■ \\x) be a finite dimensional normed linear space, (X* , \\ ■ \\x*) be its dual 
space, Z = {x\, ■ ■ ■ ,x* M } C X* \ {0}, and x* G X* . Suppose the operator x i— > (x\(x), - ■ ■ ,x* M (x)), 
x G X, is one-to-one, and the following two conditions are satisfied: (1) If x G X and x*^(x) > for 
I = 1. • • • , M , then x*(x) > 0, and (2) there exists some xq G X such that x^(xq) > for £ = 1, • • • , M . 
Then there exist nonnegative numbers Wi, £ = 1, • • • , M such that 

A I 



x*(x)=J2w e x*(x). (4.4) 
l=\ 

We will apply this proposition with X = II am for a judiciously chosen a, and x\(P) — P(xg) for 
X£ G C. The choice of a is dictated by the following Lemma H. II The inequalities of the form (|4.5p are 
known as Marcinkiewicz-Zygmund (M-Z) inequalities. 

Lemma 4.1 Let a, (3 > ~1, p > 1. There exists a constant a > depending only on a,[3,p, with 
the following property. Let C C (—1,1) be a finite p-uniform set, and m > 1 be an integer such that 
tt /(2m) < 5q < 7r/m. Then for every P G II Qm; 

ci\\P\\^a <^2\ m (fx^;z)\P(z)\ < c 2 \\P\\^ ) , 1 . (4.5) 
zee 

Moreover, if P(z) > for each z G C, then there exists c G (0, 1) (depending only on a,f3) such that 

! P(t)w a ,p(t)dt>cY / ^m(ti {aJj) ;z)P(z)>0. (4.6) 

zee 

In order to prove this lemma, we need to recall certain basic inequalities. Parts (a) and (b) can be 
found in (TTJ Theorem 4] and part (c) can be found in [211 P- 108]. We recall that we omit the mention 
of the measure p from the notations if it is the Lebesgue measure. 

Proposition 4.2 Let m > 1 be an integer, P G II m; 1 < p < oo, a, (3 > —1, and w a .p G LP . Let 

w m , a Ax) := (VT^ + l/m) 2a (VTT^ + l/m) 2/3 . (4.7) 

(a) (Markov-Bernstein inequality) 

\\P'w m ,a+i/2,0+i/2\\p < cm\\Pw a ,p\\p- (4.8) 

(b) (Nikolskii inequality) 

\\Pw m , a ^\\ r <cm^ l ' r - l '^\\Pw a ,p\\ p , l<p<r<oo. (4.9) 

(c) Let 9, ip G [0,7r], m > 1 be an integer, \9 — ip\ ~ 1/m, x = cos9, y — cosy), x < y and u,v G [x,y]. 
Then 

\ m (p {a ' p) ;u) - (l/m)w miQ+1/2:/ 3+i/2(w) ~ (l/m)w„ ha+1/2 ^ +1/2 (v). (4.10) 

PROOF of Lemma [4~T1 In this proof only, let M > 4 be an integer, C = {z k = cos9k}kLi- Since C is 
uniform, we also have qc ~ 1/m. In this proof only, let 9q — 0, 9j — 9j + qc, j : = 1, ■ • • , M — 1, 9m = tt, 
Zj = cos9j, < j < M, and Ij := [zj,Zj-\\, j — 1, • • • , M. Let n > 1 be an integer, and P G II n be 
arbitrary. We will prove first that 

jr j \ P (t) - p{ Z] )\w a , p {t)dt <cA- + (-) 2 1 yp^iii. (4.ii) 
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We observe that 1 - z 1 = 2sin 2 (#i/2) ~ l/m 2 . In view of flTO]) . 
™ Q>/3 (t)di ~ (1 - ~ ^£i!ll „ „ g^+W^ „ A m (/^U). (4.12) 



Zl „„ „. to 

Similarly, 

rZM—i 



W a ,f3(t)at ~ (1 + ZM-l) ~ 5 ~ n 



(4.13) 



In view of ()4.12|) and (14. 9|) . we have for u G 1] 



1 2 

C C CTl 

\P(u)\ I w at p(t)dt < —T\P(u)\w miat p(zi) < -^\P(u)\w m . a ^(u) < — -||P^ Q ^11 
'g m z m z 



i- 



We estimate |P(u)| J* ^ ^w a ^{f)dt (v S [—1, 5m- i]) in the same way using (|4. 13|) in place of (|4.12p . and 
deduce that 



f 1 l-ZM-i 2 

/ \P(t) - P{zi)\w a> p{t)dt + / \P(t) - P(z M )\w a ,f){t)dt < —\\Pw a Ah 

JSi J-l m 

If 2 < j < M — 1, then for t,u6 ij , we have 

w q+1/2,/3+1/2(0 ~ ^m,a+l/2,/8+l/2(*) ~ w a+l/2^+l/2(")- 

Consequently, using (|4.8p . we obtain 

M-l „ M-2 



(4.14) 



53 / |P(t) - Pfe)l«>a,/3(t)«tt < 53 / / |i y («)|du«>a 1 /s(t)«ft 
j"=l Jlj j"=2 



M-l „ „ 

^ C I3 / / li^Mha+iA/j+i/aWduCl-* 2 )" 172 * 
i=1 Ji Ji 

< cmax(0, + i - 9j) J \P\u)\w m , a+1 / 2 ,p+i/2{u)du 

CTl 

< (4.15) 

TO 

The estimates (|4.14|l and (|4.15|) together imply (|4.11|) . 

In the remainder of this proof, we consider the value of C3 fixed as in (|4.1ip . Let a = min(l, 1/(803)), 
and P £ Ham- Then with am in place of n, 



C3 

Since 



— +(— J f < c 3 (a + a 2 ) < 2c 3 a < 1/4. (4.16) 

TO \TO/ I 



/l M . M . 

\P{t)\w aJj {t)dt-Y,\ P ^)\ W *At)dt <53/ l P (*)-- P (*i)l«'a,/9( i )* 
"t j"=l "^a 7=1 



we obtain from (|4. 1 1|) that 



M 

(3/A)\\Pw a . \\i < 53 / «»a,/j(t)* < (5/A)\\Pw a ,p\\i- (4.17) 
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We now observe that for 2 < j < M — 1, and t G Ij, w a ,p(t)y/l — i 2 ~ u> w , a +i/2,/3+i/2(i)- Therefore, 
(TCTU1) implies that 

,, u , /" _ /,x Wm,a+l/2,f3+l/2(Zj) . , ( Q fl) . 

w a ,/3{t)dt~ / u; miQ , +1/2i g +1 / 2 (*)-== A m (/^ 

J/j v 1 — t m 

Together with ()4.12|) and (j4. 13[> . we have proved that 

w a ,f,{t)dt~ A m (^'«; Zj ), j = l,...,M. (4.18) 

The estimates (j4"T7|) and (|4TTg|l imply (|33]1 . 

Next, let P(zj) > for j = 1, • • • , M. In view of (|4~TTj) and (|4~T7|) , we obtain 

„i M „ 

/ FWic^pTPfe) / w a . (t)dt 

M 1 1 M /• 

< E/ l^(*)-- P (^)ka,/j(*)*< Tll^a,/3||l < / ^W*- 

.7=1 .7=1 ^ 



Consequently, 

M 



/i 2 f 

p(t)w aj p(t)dt>-y2p( Zj ) w a , 

1 ^ 7 = 1 



> 0. 



In view of (|4~18)) , this implies (|4~6j) . □ 

We are now in a position to prove Theorem 14.11 mainly using the ideas in [14j , but using a trick in 
[2"0] to prove (fC^l . 

Proof of Theorem 14.11 We let a be as in Lemma FTTI In Proposition 14. 11 we choose H am in place of 
X, \C\ in place of M, the mappings P i— > -P(z), z G C, as the set Z of linear functionals. The estimate 
(|4.5p shows that the operator P i— > (P(z)) ze c is one-to-one. We will take 

X*{P) = f P(t)d^\t) -^^m(^ h ,z)P(z), 

J- 1 zee 

where c is the constant appearing in (|4.6p . If each P(z) > 0, z G C, then (|4. 6|) implies that x*(P) > 0. 
Taking ceo in Proposition ^. II to be the polynomial identically equal to 1, we see from (|4.6j) that a;* (xq) > 0. 
Thus, all the conditions of Proposition [4J] are satisfied, and we obtain nonnegative W z , z G C, such that 



j P{t)w aiP {t) = J2( w z + cX m {^ a ' p) ;z)/2)P{z), Pen am . 

•*- x zee 

Setting w z = W z + c\ m (^ a '^; z)/2, we have proved (|4~Tj) . and also that w z > cX m (^ a '^; z)/2, z G C. 
Next, let £ G C, and n be the integer part of am/2. Then (|4.1|1 and (|2.6|) imply that 

< E ^n(M (ai/3) ;£, *) = / ^ (Q,/J) ; 6 t)i»a,/j(t) = ^(m k/5) ; e, 0- 

Hence, 

^<A n (/i^);0, £eC (4.19) 

In view of (|4.10|) and the already proved fact that ^ > (c/2)A m (/^ (Q,/3) ; £), £ G C, this implies l|4.2[) . 
The estimates (1431) are clear from (|4T2l and (14.5D. □ 



Remark. We have actually proved (|4.19|) for any positive numbers w z for which (|4.ip is valid for all 
P G IT„, whether obtained via Theorem 14. II or not. 
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5 Localized polynomial kernels 

The objective of this section is to demonstrate that in the case of the Jacobi weights, any sufficiently 
smooth, compactly supported function h, which is constant on a neighborhood of 0, is a summability 
mask of order S. This fact was essentially proved in [13], but our estimates here are more elegant. We 
prefer to generalize the summability kernels defined in Section [3] to define 

oo 

<&(/x;h,x,t) := ^ h k p k (fj,; x)p k (fi; t) , x,teR, (5.1) 

where h = {hk}^^ is a compactly supported sequence. 

The summability operator corresponding to the kernel h) is defined by 

/l oo 
$( W h, x, t)f(t)dfx{t) = hkf(M k)Pk{n; x). (5.2) 
_1 fc=0 

An interesting problem in harmonic analysis is to find conditions on the sequence h so that this operator 
is bounded in some L p ([i); i.e., ||er*(/x; h, /)|| M;P < c(h)||/|j M;p . If this is the case, the sequence h is called 
a multiplier sequence. In the important case when /i = /i^ a '^\ several conditions to ensure that h is a 
multiplier sequence are available in the literature (e.g. [THIH]). Typically, these conditions do not hold 
when p — 1 or p — oo. For our research, we therefore need to find some conditions to ensure that the 
sequence h is a multiplier sequence also in these cases. It is not difficult to see (cf. |1Q[ Proposition 2.1]) 
that this requirement is equivalent to 

sup f \$(ji,;h,x,t)\dfi(t)<c(ti). (5.3) 
set— J-i 

Theorem 5.1 Let a,f3 > —1/2, S > 1 be an integer, h k = for all sufficiently large k. Then for 
6,tpe [0,tt], 

|$(/i( a '^;h,cos0,cos^)| 

max(a,/3)+S+l/2 



^ mini (Ai + l),— — 



k=0 v ' " 



X 

_ (5.4) 

xJ2{k + l)max(a,/3)+l/2-S+m| A m /lfe |_ 
m— 1 



In particular, if h : [0, oo) — > [0, oo) is a compactly supported function that can be expressed as an S times 
iterated integral of a function of bounded variation, and h'(t) — in a neighborhood of Q, then for A > 1, 

\* x (^;h,co S e,cos<p)\ < cA^^+M^-^rnin (l, (A|g _ y|)roa L (aW+1/2 ) . (5.5) 

The following corollary is a routine consequence of Theorem 15. 11 proved in [131 Lemma 4.6] using a 
different argument. We omit the proof. 

Corollary 5.1 Let a,f3> —1/2, S > max(a,/3) + 3/2 be an integer, h k — for all sufficiently large k. 
Then 

/l S oo 

m^;h ) x ) y)\dfi^\y)<c^2^2(k + iy- 1 \^h k \. (5.6) 

- 1 3 = 1 k=Q 

In particular, if h : [0, oo) — > [0, oo) is a compactly supported function that can be expressed as an S times 
iterated integral of a function of bounded variation, then 

sup f \^ x (^ a ^;h,x,y)\dn^\y)<oc. (5.7) 
a>o, xe[-i,i] J-i 
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Analogues of (|5.5j) and (|5.7|) hold in very general situations, e.g., in the case of eigenfunctions of Laplacc- 
Beltrami operators on a smooth manifold rather than Jacobi polynomials [10] . 

Although (|5.7p is routinely proved using an estimate similar to (|5.5[) . it is sometimes possible to obtain 
such estimates without first proving a localization estimate (OH]). We have recently observed in [18] 
that there is a simple construction of exponentially localized operators from kernels satisfying analogues 
of (|5.7[) . with interesting consequences for spectral approximation of piecewise analytic functions. 

In order to prove Theorem l5.1[ we start by recalling the following estimate proved in [T31 Lemma 4.10], 
using the explicit formulas for the Christoffcl-Darboux kernel for the Jacobi polynomials: 

Proposition 5.1 Let a, (3 > —1/2, S > 1 be an integer, hk = for all sufficiently large k. Then 

|$(^);h,l,y)| 

1/2+5/2+1/4 



C < 



^minUfc + l) 2 ,^ 

S-1 

x J2(k + l) a+1 / 2 ~ m \A s ~ m h k \, ifO < y < I, 

m=Q 

oo S-1 

^(k + l)«+/ ? + 1 (* + l)~ m \A s - m h k \, if -1 < y < 0, 



v fc=0 



m=0 



(5.8) 



h,-l,y)| 

' 00 / 1 
£ mi n((fc + l) 2 ,— 



fe=0 



/3/2+S/2+1/4 



S-1 



< C ^ 



S — m j 



S-1 



j^(fc + j2 (* + ir m iA s ~ m /^ 



m=0 



i/-Ky<0, 
i/0<y<l. 



(5.9) 



Proof of Theorem [57TJ We note first that the uniqueness of orthonormalized polynomials implies that 
Pk (a ' P) (x) = {-l) k pk iM {-x), and hence, $(/i( a <0); h, x, y) = $(^' Q ); h, -x, -y). Therefore, we may 
assume that a> (3, The estimate (|5.8j) is then equivalent to (|5.4|) for the case when cos0 = 1. To extend 
the estimate for $(// a '^); h, x, y) for every x,y e [—1, 1], we recall the following product formula (|5.11[) 
by Koornwinder [8]: Let a > (3 > -1/2, K := [0, 1] x [0,tt], and for x,y E [-1, 1], r € [0, 1], u € [0,tt], let 

(1 + x)(1 + h) (1-x)(1-w) , , , - 

F(x,y;r,w) := K ^ ,yy + ^ ^ — r 2 + y/l - x 2 ^l - y 2 rcosw - 1. (5.10) 

There exists a probability measure /} = /2q,,^ on 7?. such that for n = 0, 1, • • •, and x, y G [—1, 1], 

Pn (a ' /3) (^K (a ^ ) (2/)= / ^^(iK^^^^^r,^)^^,^. (5.11) 

It follows that 

$(^;h,x,y) = [ ^ a ^-,\l,F{x, y] r,u>))djl{r,u>). (5.12) 
Jn 

Let a; = cos9,y — cost/?. Then 

Ffoyjr.w) = i(l + .x)(l + y) + -(1 - x)(l - y) + v 7 ! - i ! Vl - 

)(1 - y)(r 2 - 1) + Vl - ^V 1 - y 2 (rcosuj - 1) - 1 
= xy + Vl-^ 5 Vl-y 2 + ^ - • T )( 1 - y )( r2 _ + Vl - ^ 2 V 1 _ y 2 (rcosuj - 1), 
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and we have 



1 - F(x, y; r, u) = 1 - (xy + ^J\-x 2 ^\-y 2 ) + 1(1 - x)(l - y)(l - r 2 ) 

+ (1 — r cos bj)\J\ — x 2 \J\ — y 2 

> 1 - (xy + \/l-x 2 ^/l-y 2 ) = 1 - cos(6> - <p) = 2 sin 2 
2 



2 
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We now observe that a + S + 1/2 > f3 + 1/2, so that ijSTgjl implies 
|$0i (a ' / ' ) ;h,l, J F(i >y ;r ) w))| 

< c£ m in(( fc + l) 2 , T --i 



a/2+S/2+l/4 

X 



F(x,y;r,u) / 
J2( k + l) a+1/2 ~ m \A S ~ m h k \ 



k=0 

S-1 

X 

Q/2+S/2+1/4 



< C ^nhn((fc + 1) 2 ,^-^ 



J2(k + l) a+1 / 2 - m \A s - m h k 



k=0 N 

-> ! 

X 

m— 

Since jl is a probability measure, (|5 . 1 2[) now implies (|5.4[) . 

We observe next that if /i : [0, oo) — > [0,oo) is a compactly supported function that can be expressed 
as an S times iterated integral of a function of bounded variation, h'(t) — in a neighborhood of 0, 
and A > 1, then we may apply the above estimates with the sequence h\ whose fc-th term is given by 
h\,k = h(k/X), k — 0, 1, • • •. A repeated application of mean value theorem implies that for any s 6 K 
and integer r > 1, 

oo 

J2(k + l) s A r h x ,k < c\- s \ Arh kM < c\ s - r+1 V(h( r - 1 '>). (5.13) 

fc=0 cA<fe<ciA 

The estimate (|53|) follows from pT4"]l and (f5TT^|) . □ 



6 Some comments on computations 

There are several algorithms for computing orthogonal polynomials and the corresponding Gauss-Jacobi 
quadrature formulas, as well as expressions of the form X)fc=o a kPk(^ x) |6J. Although we have not 
necessarily used these, the localization properties of our operators and their ability to detect singularities 
and local Lipschitz exponents has been demonstrated in a number of papers, in particular, [T71 [TS1 19]. 
We offer a numerical example to illustrate the construction of quadrature formulas based on scattered 
data, using ideas described in further detail in [9J in the case of the sphere. 

Given a set C and an integer n > 1, we define a measure vc that associates the mass A„(/x;z) with 
each z e C. If {q k € Hk} is the (finite) system of orthonormal polynomials with respect to vq, there exist 
constants dj t k such that q k =^2j dj,kPj(n)- If P £ H n then with mo = \J (JbQ— 1, 1]), 



P(x)dfi(x) = P(t)J2qk(x)q k (t)dvc(tW(x) = mo P(t)J2 d oMk(t)dv c (t) 

= ^ \ m Q \ n (^; z) ^ d oMk(z) \ P{*)- 



fc=0 
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If G is the matrix with Gj t k = ^iPj{^]t)Pk{^]t)dvc{t), arid D is the matrix with (J, fc)-th entry dj t k> 
then the orthonormality of the system {%} leads to G _1 = DD T . Hence, if e = (1, 0, • • • , 0) T € R n+1 
and b = (bo, ■ • • , b n ) satsifies Gh — e, then our quadrature formula is 

/ P(x)dfj,(x) = ^2 \ m Xn(p; z)^2bjp 3 (n; z) \ P(z). 

zee [ j=o J 

This algorithm will work as long as G is a positive definite matrix. In this case, the system Gh = e can be 
solved using state of the art methods such as the conjugate gradient method. The speed of convergence 
of this method as well as the quality of the quadrature formula depends upon the condition number of 
G, which in turn, is the square of the ratio of the upper and lower bounds in the M-Z inequalities (|3 . 11 [) 
with vc in place of v and p — 2. While Gram matrices are ill conditioned in general, they seem to be 
quite well conditioned for the purpose of this application as long as the degree of the polynomial is not 
too large in relation with the number of points considered. 

As a numerical experiment, we chose 1024 points on the interval [—1, 1], chosing one point randomly 
from each of the 1024 equal subintervals of [0,7r], and taking the cosine transform to obtain points 
on [—1,1]. For various values of n, we attempted to obtain quadrature formulas exact for integrating 
polynomials in H n with respect to the Lebesgue measure, /i^ '°\ The error in the construction was 
measured by the norm of the difference between the identity matrix of order [n/2\ and the matrix of 
computed inner products Q^ z w z Pk^°'°^ {z)Pe^°'°^ ( z )) , k, £ = 0, ■ • ■ , [n/2\ . Each computation was carried 
out 30 times using Matlab 7.1 running on Intel Core 2 CPU with 2.13GHz, and 3GB RAM. The averages 
of the various parameters of interest are summarized in Table [TJ It is interesting to note that we were able 
to obtain suprizingly accurate quadrature formulas even for degree 1023. Even though there were several 
negative weights in this case, the total variation of the resulting discrete quadrature measure seems to 
be surprisingly small. 



n 


cond 


pos 


EM 


error 


time 


256 


1.5342 


1024 


2.0000 


2e- 14 


0.4769 


512 


3.3536 


1023.92 


2.0000 


9e- 14 


1.2441 


768 


25.2965 


1012.28 


2.0088 


2.4e- 13 


2.6159 


896 


76.8983 


945.60 


2.1471 


3.8e- 13 


3.7419 


1023 


11270.5423 


885.64 


2.6715 


2.75e- 12 


9.5317 



Table 1: Averages for 30 trials of: n — degree of exactness, cond is the condition number of the original 
Gram matrix, pos is the number of positive weights, error is the norm of the difference between the 
matrices as described in the text, (xe-m means x* 10 _m ), time is the time required for each computation, 
in seconds. 
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